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Engaging strongly resonant interactions allows dramatic enhancement of functionalities of many 
electromagnetic devices. However, resonances can be dampened by Joule and radiation losses. While in 
many cases Joule losses may be minimized by the choice of constituting materials, controlling radiation 
losses is often a bigger problem. Recent solutions include the use of coupled radiant and sub-radiant modes 
yielding narrow asymmetric Fano resonances in a wide range of systems, from defect states in photonic 
crystals and optical waveguides with mesoscopic ring resonators to nanoscale plasmonic and metamaterial 
systems exhibiting interference effects akin to electromagnetically-induced transparency. Here we 
demonstrate theoretically and confirm experimentally a new mechanism of resonant electromagnetic 
transparency, which yields very narrow isolated symmetric Lorentzian transmission lines in toroidal 
metamaterials. It exploits the long sought non-trivial non-radiating charge- current excitation based on 
interfering electric and toroidal dipoles that was first proposed by Afanasiev and Stepanovsky in [J. Phys. A 
Math. Gen. 28, 4565 (1995)]. 



Toroidal dipole (and higher toroidal multipoles) is not a part of the standard multipole expansion that is 
usually presented as a sum of contributions from two families of conventional elementary radiation sources, 
namely the magnetic and electric dynamic multipoles 1 ' 2 . In a spherical coordinate system, magnetic multi- 
poles are defined by transversal components of oscillating current density (e.g. currents flowing on a surface of a 
sphere), while electric multipoles are attributed to oscillating charge density. Fields of an arbitrary charge-current 
configuration also contain a contribution from oscillating radial components of the current density, which is, 
however, dropped in the commonly used long- wavelength approximation rendering the standard multipole set 
incomplete, as was pointed out in Refs. 3-5. These radial components of the current density give rise to an 
additional independent family of elementary radiation sources, the dynamic toroidal multipoles, which are 
physically rather than just semantically different from the conventional magnetic and electric multipoles 6 ' 7 . 
While electric and magnetic dipoles can be visualised as a pair of opposite charges and a circular current 
respectively, toroidal dipole results from currents j flowing on the surface of a torus along its meridians (poloidal 
currents), as illustrated in Fig. la. Its moment is given by 



1 

10c 



[(r-j)r-2r 2 j] d\ 



where r is the coordinate vector with its origin placed in the center of the torus 3 . 

In physics of nuclei static toroidal moment is a well-established notion, also referred to as the anapole. It has 
been calculated for a number of nuclear systems and investigated in details theoretically and experimentally, in 
particular, for its role in nuclear-spin-dependent atomic parity violation (see 8 9 and references therein). Toroidal 
moment has also been identified from first principle calculations in carbon macromolecules 10 and ferroelectric 
nanostructures 11 . However, the dynamic excitations of toroidal dipole and higher toroidal multipoles were 
overlooked for some time in classical electrodynamics as their manifestations are often masked by much stronger 
electric and magnetic multipoles. Only recently, a spectrally isolated transmission resonance corresponding to the 
toroidal excitation has been observed in an artificially engineered electromagnetic material (metamaterial), where 
contributions from electric and magnetic multipoles were deliberately suppressed by design 12 . Since then micro- 
wave and optical toroidal resonances were identified in a number of metamaterial and plasmonic systems 13 " 16 . 
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Figure 1 | Toroidal metamaterial. (a) Poloidal currents flowing on a surface of a torus along its meridians create toroidal dipole moment T. (b) Metal 
screen with a dumbbell- shaped aperture is the structural element of toroidal metamaterial. Dashed arrow m represents axis of its mirror symmetry, 
(c) Photograph of the assembled metamaterial slab, an array of 15 X 16 toroidal aperture-based metamolecules. Inset shows a close-up view of one of the 
array's column with 8 -fold symmetry. 



Importantly, the emissions of oscillating electric and toroidal 
dipoles have the same angular distribution and parity properties, 
although there are fundamental physical differences between the 
two elementary excitations: (i) radiated power for electric and toroidal 
dipoles scales as correspondingly 5 " 7 w 4 and w 6 ; (ii) electric dipole is 
induced by the electric field of the wave, E, while toroidal dipole is 
induced by 12 curl B ~ dE/dt; (iii) vector potential fields A generated 
by oscillating electric and toroidal dipoles are different and the dif- 
ference cannot be removed by gauge fixing 17 . Thus, when the poloidal 
current mode (toroidal dipole) and bipolar charge distribution (elec- 
tric dipole) are collocated and coherently oscillate at the same fre- 
quency, their radiated electromagnetic fields can be made to interfere 
destructively and disappear everywhere outside the resulting charge- 
current configuration, as was first pointed out by Afanasiev and 
Stepanovsky 17 and demonstrated with further numerical studies 18 . 
This gives an interesting opportunity of creating artificial electromag- 
netic materials, which exhibit a very narrow isolated transparency 
band when the radiation losses are suppressed through the interfer- 
ence of toroidal and electric dipolar emission. The underpinning 
mechanism of the transparency and resonant dispersion here are 
somewhat different from those typically encountered in the plasmo- 
nic and metamaterial systems mimicking electromagnetically- 
induced transparency, where the interference between two spectrally 
separated multipolar modes produces a sharp asymmetric Fano res- 
onance on the background of another much wider resonant band 19 . 

Results 

We demonstrate this novel mechanism of radiation suppression for 
metamaterials based on a dumbbell-shaped aperture element. It is 
made in a thin metal plate and resembles a meridianal cross-section 
of a toroidal void (see Fig. lb). This is a very special electromagnetic 
system in which the incident wave E(f) = E 0 e ioJt polarized parallel to 
the symmetry axis of the aperture m induces both electric P and 
toroidal T dipolar moments. Indeed, the electric field of the wave 
drives charge separation p (t) = p 0 e icot across the waist of the aper- 
ture, which gives rise to an oscillating electric dipole moment 



P(t) = Ye ia)t ~ mp 0 e ioJt oriented along m (see Fig. lb). The charge 
displacement also gives rise to counter-rotating (poloidal-like) cur- 
rents j(t) oscillating along the edges of the circular cuts, as shown 
in Fig. lb. The poloidal currents are proportional to the time deriv- 
ative of the charge displacement and produce toroidal moment 
T(t) = Te ia>t ~ miwp Q e l0Jt that is also oriented along m and oscil- 
lates coherently with the electric dipole lagging a quarter of the 
period behind the latter. Now, if the field radiated by the electric 
dipole E P (t) is proportional to the second time derivative of P(t), 
E P (t) ~ icopo e ioJt , the field radiated by the toroidal dipole is propor- 
tional to the third time derivative of 17 T(t), i.e. Ej{t) ico 3 Po eicof - 

Therefore E P (t) and E T (t) scattered by the dumbbell- shaped aperture 
are in anti-phase and always interfere destructively (see Supplemen- 
tary Materials for more details). The amplitude of the poloidal cur- 
rents increases resonantly when the wavelength of incident radiation 
X = 2n c/co becomes close to the circumference of the aperture — 4tt R 
(where R is the radius of the circular cuts) leading to the enhance- 
ment of toroidal emission, which in principle can be made to cancel 
electric dipolar scattering completely. 

In reality, the incident electromagnetic wave induces in such a 
structure higher oscillating multipoles as well, most notably mag- 
netic quadrupole Q (m) , which also results from the pair of counter- 
rotating currents. The scattering contribution from this multipolar 
current mode can be effectively suppressed in an aperture-based 
structure of higher rotational symmetry such as, for example, 4-fold 
or 8-fold symmetric toroidal metamolecules shown in Figs. 2a and 
2b. We demonstrated this by modeling numerically the interaction of 
the metamolecules assembled in two-dimensional arrays (slabs of 
toroidal metamaterials) with a normally incident linearly polarized 
plane wave. Calculated densities of the induced electrical currents 
were used to compute scattered powers of the conventional multi- 
poles and toroidal dipole associated with each metamolecule (see 
Supplementary Materials for details). 

Our calculations showed that, for example, in the metamaterial 
with 4-fold symmetric metamolecules the emission of the standard 
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Figure 2 | Multipole excitations in toroidal metamaterials. Panels (a) and (b) show metamaterial's unit cell with 4-fold (a) and 8-fold (c) symmetry; 
a = 16.5 mm, R = 2.5 mm, w = 1.2 mm, separation between the centers of the circular cuts is 7.0 mm. Panels (c) and (d) show dispersions of multipolar 
scattering rates calculated for 4 strongest multipoles induced in 4-fold (c) and 8-fold (d) symmetric metamolecules. Arrows indicate locations of the 
corresponding transparency resonances. 



multipoles other than electric dipole P is small: close to the resonance 
(which is at around 1 1 GHz for a centimetre-sized metamolecule) 
electric quadrupole Q (e) and magnetic dipole M exhibit scattering 
rates that are factors of 10 and 30 smaller than that of P, while 
Q (m) is a factor of 10 5 weaker here (see Fig. lc). Toroidal dipole 
provides the second strongest contribution at the resonance being 
responsible for more than 30% of the total scattering. Its presence can 
be detected in the near-field as closed loops of magnetic field-lines 
confined within the metamolecule and threading the circular sec- 
tions of the apertures (see left inset to Fig. 3). 

Our calculations also show that at the resonance the complete 
destructive interference of radiated fields take place in a loss-less 
toroidal metamaterial (metals can be treated loss-less in the micro- 
wave part of the spectrum). This corresponds to the total transpar- 
ency of the toroidal metamaterial. For instance, for the 4-fold 
symmetric system the prevailing electric dipolar scattering is can- 
celled on about 65% by toroidal scattering and on 35% by scattering 
from other multipoles (see Supplementary Materials for more 
details) creating a background- free resonance transparency peak with 
symmetric Lorenzian profile and Q = 35, which is already quite high 
by the standard of the microwave metamaterials (see Fig. 3). In the 
metamaterial composed of 8 -fold symmetric metamolecules electric 
and toroidal dipolar scattering is mutually cancelled on 96%, leading 
to a transparency peak with Q ~ 520. The increased role of toroidal 
dipolar mode in the 8 -fold symmetric structure is also evident from 
the calculated magnetic near-field map (compare the insets to Fig. 3). 

We confirmed the appearance of the transparency resonances due 
to electric-toroidal dipolar interference with experiments conducted 



in the microwave part of the spectrum. For that we constructed 
toroidal metamaterials from strips of thin stainless steel plates with 
dumbbell-shaped apertures, which were assembled into 4- or 8-fold 
structures (see Fig. lc). The experimental data shows a very good 
qualitative agreement with the theory reproducing fully the main 
features of the metamaterials' response predicted by the modelling: 
the appearance of narrow isolated bands of transparency with sym- 
metric Lorentzian profile and very large Q- factor, which reaches a 
record-high value of 320 in the case of 8-fold symmetric metamater- 
ial (see Fig. 3). The slight shift of the observed resonances towards 
higher frequencies is attributed to the shortening of the outer effec- 
tive circumference of the apertures due to the non-zero thickness of 
the metal strips. For the strip thickness of 0.8 mm our estimate 
gives the relative shortening and correspondingly blue shifts of about 
5% and 14% for 4- and 8-fold symmetric metamolecules respectively, 
which agrees well with the results of the actual measurements 
(see Supplementary Materials for details). Somewhat lower Q- 
factors measured for the toroidal metamaterials result from the 
inhomogeneous broadening of the resonances due to fabrication 
tolerances and imperfections of the sample assemblies. The latter is 
also responsible for the incomplete transparency, as well as the 
appearance of the broad asymmetric pedestal of the high-fidelity 
resonance. For example, given 0.2% variation of the resonance fre- 
quency among 8 -fold symmetric metamolecules in the assembled 
array (which correspond to the fabrication tolerance of 0.06 mm), 
the amplitude of its transmission resonance would reduce to 4% with 
the Q-factor dropping to about 300. Other factors that limited the 
level of transparency in the experiments are residual curvature of the 
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Figure 3 | Transmission response of toroidal metamaterials. (a) Blue and green solid curves show simulated spectra of metamaterial slabs based on 4- 
fold and 8-fold symmetric metamolecules respectively, dots - experimentally measured data (multiplied by the factor of 50 in the case of 8-fold symmetric 
structure). Orange-filled areas represent symmetric Lorentz line profiles fitting simulation and experimental data for the 4-fold symmetric structure. 
Insets show simulated distribution of magnetic field lines and magnetic field intensity I H I corresponding to resonantly excited modes in 4-fold (left) and 
8-fold (right) symmetric metamolecules. (b) Simulated (solid line) and measured (dots) spectra of 8-fold symmetric metamaterial fitted by symmetric 
Lorentz line profiles (orange-filled areas). 



incident wavefront and diffusive scattering from the metamaterial 
structures. 

Discussion 

We point out that simultaneous presence of collocated electric and 
toroidal dipolar excitations in our metamaterials may create a unique 
situation, where vector potential A generated outside the arrays is 
non-zero even though the re-emitted (in particular, reflected) elec- 
tromagnetic fields are completely cancelled via the destructive 
interference. Indeed, although the electromagnetic field emission 
characteristics of oscillating P(f) and T(f) are identical, the vector- 
potential fields they produce are quite different 17,20 . Under the con- 
dition of total coherent cancellation of their radiation the resulting 
vector potential propagating in the far-field is given by 

A(r,t) = - r -^(rT)e- ikr+i0Jt . 

The complete expression for A featuring also the near-field con- 
tributions can be found in Supplementary Materials. 

The discussion on the independent physical significance of the 
vector potential started from the famous paper by Aharonov and 
Bohm 21 and has developed into an active research topic that attracted 
a large number of publications (see, for example 22-24 and references 
therein). In particular, it was theoretically shown that the Aharonov- 
Bohm effect due to time-varying vector potential of a coherent light 
source could be used for quantum-nondemolition detection of 
photons emitted by the source 25 . Intriguingly, Afanasiev and 
Stepanovsky suggested in their paper 17 that time- dependent electro- 
magnetic potentials produced by their non-trivial non-radiating 
system of collocated electric and toroidal dipolar excitations "... 
can be used as a new channel for information transfer (by modulat- 
ing the phase of the charged-partial wavefunction) and for the per- 
formance of time -dependent Aharonov- Bohm-like experiments.". 
Those predictions were disputed by Marengo and Ziolkowski, who 
argued "... that electrodynamically not only the fields but also the 
associated potentials are unobservable everywhere in the exterior of a 
spatially localized non-radiating source" and "... Aharonov- Bohm 
effects associated with non-radiating potentials are possible only in 
static situations." 26 . Resolving this apparent controversy and further 
understanding the role of time- dependent electromagnetic potentials 



would require practical realization of non-trivial non-radiating 
sources, which have so far remained merely a peculiar theoretical 
concept. 

We believe that the observed resonant transparency phenomenon 
is the first example of the manifestation of the non -trivial non- radi- 
ating charge-current excitation, while toroidal metamaterial arrays 
operating at the transparency resonance could serve as the source of 
the non-compensated vector-potential. Even though in our case the 
destructive interference between electric and toroidal dipole 
moments does not lead to the complete cancellation of their scatter- 
ing, the residual dipolar radiation (which converges in forward and 
backward directions for an infinite array of sub -wavelength peri- 
odicity) is suppressed by higher multipoles, while the gauge-invari- 
ant vector-potential should be generated by the mutually 
compensated contributions of P and T. Here, due to the two-dimen- 
sional periodicity of the metamaterial arrays, the non-trivial com- 
ponent of the vector-potential field is expected to localize in the plane 
of the metamaterial structure. 

In conclusion, we have identified a class of metamaterials support- 
ing a novel mechanism of resonant electromagnetic transparency, 
which is different from Fano interference typically encountered in 
the plasmonic and metamaterial systems mimicking electromagnet- 
ically-induced transparency. It exploits destructive interference 
between spatially and spectrally collocated and coherently oscillating 
induced electric and toroidal dipoles, and produces very narrow and 
characteristically symmetric Lorentzian transparency lines with Q- 
factors exceeding 300. Such mode of the metamaterial resonant 
excitation corresponds to a long-awaited implementation of the 
non-trivial non-radiating charge- current configuration that might 
generate waves of gauge-irreducible vector potential in the absence 
of scattered (reflected) electromagnetic fields. 

Methods 

Numerical modelling. The metamaterial slabs were described through periodic 
boundary conditions applied to the corresponding unit cell's facets in X and Y 
directions, as indicated in Figs. 2a and 2b. The dumbbell- shaped aperture resonators 
forming 4-fold and 8-fold symmetric toroidal metamolecules were assumed to be cut 
in infinitely thin sheets of perfect electric conductor in accordance with the design 
specifications. Electromagnetic response of the metamaterial slabs was simulated 
using commercial full three-dimensional Maxwell equations solver based on the finite 
element method, COMSOL 3.5a. The simulations also provided data on the densities 
of electrical currents induced in the metamaterials by the incident wave, which was 
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used to compute scattered powers of the conventional multipoles and toroidal dipole 
associated with each metamolecule. 

Samples fabrication. Metamaterial slabs were constructed from strips of thin 
stainless steel plates with dumbbell- shaped apertures. The stainless steel strips, which 
had the thickness of 0.8 mm, were patterned using chemical etching and assembled 
into columns of 4 -fold and 8 -fold symmetric structures. The metamaterial slabs were 
formed by 15 such columns, which contained 16 toroidal metamolecules each (see 
Fig. lc). All dimensions of the design features were identical to those used in our 
modelling (Fig. 2). 

Microwave measurements. Transmission spectra of the metamaterial samples were 
measured in the Emerson & Cuming microwave anechoic chamber using vector 
network analyser (Agilent E8364B) and a pair of broadband linearly polarized horn 
antennas (Schwarzbeck Mess-Elektronik BBHA 9120D) equipped with dielectric lens 
collimators. 
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